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Abstract 

The dynamical properties of a 2D Heisenberg model with dipolar interactions and perpendicular anisotropy 
are studied using Monte Carlo simulations in two different ordered regions of the equilibrium phase diagram. 
We find a temperature defining a dynamical transition below which the relaxation suddenly slows down and the 
system aparts from the typical Arrhenius relaxation to a Vogel-Fulcher-Tamann law. This anomalous behavior is 
observed in the scaling of the magnetic relaxation and may eventually lead to a freezing of the system. Through 
the analysis of the domain structures we explain this behavior in terms of the domains dynamics. Moreover, we 
calculate the energy barriers distribution obtained from the data of the magnetic viscosity. Its shape supports our 
comprehension of both, the Vogel-Fulcher-Tamann dynamical slowing down and the freezing mechanism. 
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1. Introduction 

Thin magnetic films have been the subject of 
intense attention over the last two decades [1-3]. 
These studies have been motivated both by the 
technological importance of these systems and by 
the insight that they provide into the fundamental 
role of interactions at the atomic level. Advances in 
film growth techniques and experimental control of 
spin-spin interactions have impulsed this interest. 

It is well established now the fundamental role 
of competing interactions in the emerging fea- 



tures of low dimensional systems. Among a wide 
scope of theoretical and numerical investigations 
on equilibrium [4-6] and dynamical properties [7-9] 
of several model Hamiltonian of low-dimensional 
magnetic systems, Heisenberg-like models are one 
of the most intensively used to approach real mag- 
netic materials. [10,4] In fact, anisotropic versions 
of the Heisenberg model describe several com- 
pounds as BaCo 2 (As0 4 ) 2 ,[ll] CoCl 2 -GIC[12] and 
Rb 2 CrCl 4 [13] 

Some effort in the study of the non-equilibrium 
regime of closely related models has been done[14- 
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17] and a detailed characterization of the equilib- 
rium features of the Heisenberg Hamiltonian con- 
sidering a dipolar long-range repulsive term and 
a varying perpendicular anisotropy were depicted 
by numerical and analytical methods in the last 
few years. [4, 18, 19] In particular, the transition in 
which spins reorient from an in-plane ferromag- 
netic order to an out-plane perpendicular stripes 
with decreasing temperature, already seen in ex- 
perimental works, [20] was encountered in a well- 
localized anisotropy range. However, a complete 
characterization of this system is far from been es- 
tablished, and that is particularly true in relation 
with its dynamical behavior. In this paper we want 
to make a step forward in this direction. 

Through Monte Carlo simulations we study the 
dynamics of the model analyzing the magnetic re- 
laxation behavior for several values of tempera- 
ture. We explore dynamical properties in two dif- 
ferent ordered regions of the phase space near the 
Spin Reorientation Transition. We show that there 
is a temperature T* below which the dynamic slows 
down. This behavior is present for the two phases 
around the reorientation transition. We explain 
this behavior in terms of non-Arrhenius relaxation 
due to complex domain dynamics. 

The paper is organized as follows: In section 2 
we present the model making a review of some of 
its more relevant properties and give details about 
the Monte Carlo simulations. Then, in section 3 
we present and discuss our results understanding 
the anomalous behavior of magnetic relaxation by 
means of the analysis of the domains dynamics and 
the energy barriers associated to the different re- 
gions. Finally, section 4 is devoted to summarize 
our results and present the conclusions of the work. 



2. Model and simulation 

We consider a two-dimensional square lattice of 
Heisenberg spins of unit modulus \si\ — 1 with an 
anisotropy perpendicular to the plane of the lat- 
tice, interacting through the dimensionless Hamil- 
tonian 

1 (ij) 



ijtj \ *J ij J 

where the exchange S and anisotropy constants rj 
are normalized relative to the dipolar coupling con- 
stant, <i,j> stands for a sum over nearest neigh- 
bors pairs of sites in the lattice, (i, j) stands for a 
sum over all distinct pairs and tv, = r*i — rj \ is the 
distance between spins i and j. 

Recently, Carubelli et all [4] determined the r\-T 
phase diagram of this model by means of Monte 
Carlo simulations. Maintaining the same relative 
strength between exchange and dipolar interac- 
tions they showed the existence of two different 
ordered phases at low temperature. A planar fer- 
romagnet or a perpendicular striped phase estab- 
lishes depending on the value of rj. Moreover, the 
diagram obtained in [4] qualitatively reproduces 
detailed measurements of the reorientation tran- 
sition encountered experimentally in samples of 
Fe/Ni/Cu(001)[20] suggesting that the film thick- 
ness acts as an effective inverse anisotropy. 
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Fig. 1. Scheme of the phase diagram of a system of size 
32 X 32 for 5 = 3. The dynamical calculations are done near 
the Spin Reorientation Transition (shown as a bold line). 

We focused on the study of quenches for a sys- 
tem of size N — 48 x 48 and S = 3 from full out-of- 
plane magnetized configurations into planar and 
striped ordered regions. Several temperatures be- 
low T — 0.9 corresponding to the ordered phases 
were explored for two anisotropy values 77 = 5.5 
and 7] = 7.0. In figure 1 a scheme of the phase dia- 
gram is shown and the Spin Reorientation Transi- 
tion is marked with a bold line. 
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All the simulations were done using the 
Metropolis algorithm and infinite periodic bound- 
ary conditions were imposed on the lattice by 
means of the Ewald sums technique. The trial for 
the spin flip was chosen random orientation 
uniformly distributed in the sphere. 

The time evolution of the out-of-plane mag- 
netization M z = 1 /N J2 i sf was recorded start- 
ing from the saturated configuration s\ = 1. 
The behavior of this observable is in all cases a 
monotonous decay to zero due to the lack of total 
magnetization of the striped configurations as well 
as the complete in-plane orientation of spins in the 
planar-ferromagnetic phase. Typical time orders 
in Monte Carlo steps were t = 10 m.c.s. and av- 
erages over thermal noise involved a few hundred 
of realizations. In the quench to the striped phase 
rj = 7.0, it was impossible to relax the system at 
temperatures below T = 0.4. There was not clear 
magnetization decay up to t — 10 7 m.c.s. 



3. Results and Discussion 

We start showing in figure 2 the out-of-plane 
magnetization relaxation corresponding to both 
quenches at different temperatures. 

As can be seen there is a central difference be- 
tween quenches in the two ordered regions. For 
i] = 7.0 the ground state of the system is a striped 
configuration and, as we already mentioned, there 
exists a temperature T — 0.4 below which the dy- 
namics completely freezes. Contrary to the 77 = 5.5 
case where degeneration of planar ferromagnetic 
ground state facilitates domain relaxation, for rj = 
7.0 the dynamical flipping of spins corresponding 
to adjacent striped domains is naturally slower. 

The scaling of M z for 77 = 5.5 at high tempera- 
tures is shown in figure 3 as a function of T\n(t/r n ). 
There was no possible scaling in this functionality 
for lower temperatures. 

It means that below T < T* — 0.5 the system 
does not longer relax following the Arrhenius law. 
[21] In fact, despite many works validating the ap- 
plicability of the scaling of figure 3 to interacting 
systems, [22,16] there is no guarantee that inter- 
acting spins may be always described in terms of 
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Fig. 2. Out-plane magnetization for several temperatures 
with rj = 5.5 (up) and rj = 7.0 (down). 
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Fig. 3. Scaling of magnetization for several temperatures 
with r\ = 5.5. 

thermal activation of the Arrhenius type. 

The most simple and common extension for 
the Arrhenius equation is the well known Vogel- 
Fulcher-Tammann (VFT) law which is given by 



t(E) = T e T - T o 



(2) 
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where the new parameter To is known as the VFT 
temperature. This functionality provides usually 
a good description of experimental and numerical 
data and is the salient phenomenological descrip- 
tion for the a-relaxation dynamics of glass-forming 
systems. [23] 

Although the microscopic origin of the VFT law 
is still under debate, it is generally accepted that 
a clue feature of it is the notion of cooperativity in 
the movement of the elementary constituents of the 
system. Thus, many theories of the glass transition 
are based on the assumption that cooperatively re- 
arranging regions (domains or clusters) merges at 
low temperatures. Inside them magnetic moments 
move in some concerted way that distinguishes 
these regions from their surroundings. [24,25] 
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Fig. 4. Scaling of magnetization for all temperatures at 
rj = 5.5 considering a VFT relaxation law. Collapse is 
obtained for tq = 7 and T> = —0.5. 

From expression (2) it is direct to see that the 
expected scaling variable now takes the form (T — 
T )ln(t/ro). In figure 4 this scaling is shown for 
temperatures below T = .5. On the contrary, there 
is no collapse for higher temperatures where the 
data follows the typical Arrhenius relaxation seen 
in figure 3. It is worth to note that the value of 
the parameters for the best collapse are shown in 
the figure and involve a negative VFT temperature 
T = —0.5. As far as we know, this feature has 
being observed only in a few materials[26,27] and 
its interpretation remains unclear. 

In order to confirm this unusual feature we stud- 
ied the relaxation times t(T) defined as the time 
at which the out-plane magnetization decays to a 
value M z = 0.01. In figure 5 we see the Arrhenius 
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Fig. 5. Relaxation times of the sample magnetization for 
temperatures inside the non- Arrhenius regime. The solid 
line corresponds to the best fit with the VFT law (2). 
Parameter values are to = 35 ± 8, To = —0.47 ± 0.03 and 
U = 3.3 ±0.3. 

plot of t and its fit using the VFT form. As can be 
noticed there is a very good agreement between the 
data extracted from the simulation and the expres- 
sion (2). Confirming our previous results, the VFT 
temperature extracted from this fitting is To = 
-0.47 ±0.03. 
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Fig. 6. Scaling of magnetization for several temperatures 
with r) = 7.0. 

For the r\ = 7.0 case figure 6 shows that the 
Arrhenius-like relaxation holds down to a temper- 
ature T* = 0.4. The data corresponding to T = 
0.4 was impossible to collapse in the same figure 
and, since it is the only relaxation data below the 
Arrhenius zone, a proper collapse considering ex- 
pression (2) can not be found. The freezing of the 
system is stronger that for the case r\ — 5.5. 
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3.1. Cluster Analysis 

Now it is worth to do a more deep study to un- 
derstand the nature of the slow-down encountered 
that makes the system to apart from the Arrhe- 
nius law and may eventually freeze the relaxation 
dynamics. Firstly, the disagreement with the loga- 
rithmic scaling and better adjusting to a VFT re- 
lation reflects an intrinsic physical phenomenology 
of the model under study. That is, the emergence 
of larger times for the cooperative relaxations as 
temperature is decreased below T* . Moreover, the 
dynamical freezing for rj — 7.0 could be a feature 
coming from the combined action of the slow non- 
Arrhenius relaxation and the context of a ground 
state of perpendicular-oriented spins that natu- 
rally develops higher relaxation times due to higher 
energy barriers between local minima. 

Due to the three-dimensional nature of the spins 
in the system, the visual inspection of the config- 
urations involves two type of formats: the in-plane 
and the out-plane projection. For the in-plane pro- 
jection of the spins we use the format shown in the 
left side of figure 7. There, the spin is represented 
with a single arrow and the dark intensity is pro- 
portional to the modulus of the projection, so that 
white corresponds to a completely out-plane spin. 
On the other hand, for the out-plane projection 
the format is seen in the right side of the same fig- 
ure. In this case every spin is a tiny square and the 
modulus of the out-plane projection goes from — 1 
(blue) to 1 (red), so that white corresponds to a 
completely in-plane spin. 

Figure 7 shows a typical evolution of the sys- 
tem with r\ = 5.5 and T = 0.6. Starting from the 
fully out-plane magnetized state, the spins begin 
to depart from this position, some of them even 
inverting its original direction. At some point sev- 
eral ground-state in-plane domains merges, these 
domains are eventually separated by out-plane 
structures that continuously disappears as the 
time grows and equals r and the system orders as 
a planar ferromagnet. 

In order to compare this evolution for several 
temperatures we have to be careful. Since even in 
the Arrhenius regime (T > T* ) it is clear that the 
relaxation slows with decreasing temperature, di- 




Fig. 7. Typical system configurations for the r\ = 5.5 quench 
at T = 0.6 and several times. The left (right) side shows 
the in-plane (out-plane) projection of the spins. From up 
to down correspondent times are t = 64, 256, 1024, 16385. 

rect inspections to complete evolutions above and 
below T* does not leads to useful conclusions. In- 
stead, what we do is to fix certain equivalent stage, 
M z = 0.5, in the relaxation process at different 
temperatures and then inspect the typical config- 
urations. 

Figures 8 and 9 are typical configurations above 
and below T* respectively for which M z — 0.5, 
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Fig. 8. Typical system configurations for the rj = 5.5 quench 
at M z = 0.5 and several temperatures above T*. From up 
to down T = 0.8, 0.7, 0.6, 0.5. 

that is, in an equivalent macroscopic stage of the 
relaxation. The figures show several temperatures 
and again it is clear that the type of relaxation 
is different between the high and low temperature 
regions. Above T* the out-plane magnetization is 
mainly supported by uncorrelated spins that re- 
mains with a non-zero out-plane projection. On 
the other hand, below T* the same magnetization 




Fig. 9. Typical system configurations for the r\ = 5.5 quench 
at M z = 0.5 and several temperatures below T*. From up 
to down T = 0.4, 0.3, 0.2, 0.1. 

is due to well-structured out-plane reminiscences 
that exactly corresponds with the defects of the 
cooperative in-plane ferromagnetic clusters. That 
is, for lower temperatures the depart from out- 
plane saturation is characterized by the movement 
of highly correlated domains of in-plane spins or, 
what is the same, by the slow movement of well- 
defined (out-plane) defects of these in-plane struc- 
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Fig. 10. Typical system out-plane spin projections ranging 
from —1 (black) to 1 (white) for the r\ = 7.0 quench at 
T = 0.8 (left) and T = 0.4 (right). From up to down 
correspondent times are t = 64, 256, 1024, 16385. 



to reach the ground state. For higher temperatures 
the system is able to relax to a configuration in 
which the correct striped domains merges. On the 
other hand for T — 0.4 the evolution is frustrated 
by an abnormal growth of the stripes width. 

3.2. Energy barriers distributions 

In previous works [22, 16] it has been shown that, 
in the range of validity of the Tln(t/ro) scaling, the 
effective distribution of energy barriers (/&(£?)) 
contributing to the long time relaxation process 
can be obtained from the magnetic relaxation 
collapse simply by performing its logarithmic 
time derivative dM (t) / d\n(t) . That is, studying 
the magnetic viscosity S. This distribution fb{E) 
represents a time independent distribution that 
would give rise to an identical relaxation curve. 
This quantity is calculated in this subsection from 
the numerical data. 

Moreover, non-Arrhenius relaxation following 
equation (2) can be seen as a simple rescaling of 
T. Thus, since the derivation of S ~ fb(E) does 
not involve temperature derivatives, [28] the terms 
remain unchanged and the relation between mag- 
netic viscosity and energy barriers distribution 
still holds. So we are able to obtain the statisti- 
cal information of the energy barriers above and 
below the dynamical transition temperature. 
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Figure 10 shows typical out-plane configurations 
for the quench at rj — 7.0. In this case in-plane 
views shows no order and no significant spin pro- 
jections and therefore are uninteresting. As can be 
seen, the structure of the out-plane projections at 
T = 0.4 are consistent with an increasing difficulty 



Fig. 11. Energy barriers distribution explored for the sys- 
tem with 7] = 5.5 above and below T* . 

The results shown in figure 11 for the two dif- 
ferent dynamical regimes at r\ — 5.5 confirms our 
expectations. Indeed, the slower non-Arrhenius re- 
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laxation involve higher energy barriers due to the 
cooperative nature of the domains movement. 

In the 77 = 7.0 case some light can be spread 
on its anomalous behavior. If the energy barriers 
the system has to overcome are high enough, a 
Monte Carlo out-equilibrium dynamics freezes in 
some metastable configuration. So, it is worth to 
comparing the form of ft(E) between the quench 
to planar ferromagnetic phase and to the perpen- 
dicular striped region. Figure 12 shows the result 
of analysis of energy barriers in the two regions in 
which the Arrhenius relation holds. 
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Fig. 12. Energy barriers distribution explored for the sys- 
tem at T > T* at both sides of the Spin Reorientation 
Transition. 

As can be seen, the distribution associated to 
the striped phase involves higher values of energy 
barriers than that of the planar ground state. The 
shape of the distribution in the n — 7.0 case is also 
wider suggesting a heterogeneous distribution of 
domains. 

This picture supports our believe that, once the 
system has formed striped structures, the energy 
barriers to be overcome are higher. Indeed they are 
so high that they can eventually freeze below T*. 



4. Conclusions 

Through numerical Monte Carlo simulations we 
have studied the dynamical properties of a two di- 
mensional Heisenberg model with dipolar interac- 
tions and perpendicular anisotropy. In particular, 
the analysis of the out-of-plane magnetization dy- 



namics suggests the existence of a dynamical tem- 
perature transition T* . Below this temperature the 
dynamics changes from an Arrhenius to a VFT 
type (77 = 5.5) or the system freezes (rj = 7.0). 

The interpretation of this anomalous slow down 
is understood as a result of the change from the 
single-spin elementary movement to a new mech- 
anism based on the coherent movement of the ori- 
ented ground-state domains. The same explana- 
tion can be done in terms of the pinning of the de- 
fects merging in the frontiers of ground-state do- 
mains of different orientations. This explanation is 
confirmed by the analysis of typical configurations 
during the dynamics and by the shape of the en- 
ergy barriers distributions obtained studying the 
magnetic viscosity of the model above and below 
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